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ABSTRACT 



Aims. The properties of the early-type stars in the core of the Westerlund 2 cluster are examined in order to establish a link between 
the cluster and the very massive Wolf-Rayet binary WR20a as well as the Hn complex RCW49. 

Methods. Photometric monitoring as well as spectroscopic observations of Westerlund 2 are used to search for light variability and to 
establish the spectral types of the early-type stars in the cluster core. 

Results. The first light curves of the eclipsing binary WR 20a in B and V filters are analysed and a distance of 8 kpc is inferred. Three 
additional eclipsing binaries, which are probable late O or early B-type cluster members, are discovered, but none of the known early 
O-type stars in the cluster displays significant photometric variability above 1% at the 1-tr level. The twelve brightest O-type stars are 
found to have spectral types between 03 and 06.5, significantly earlier than previously thought. 

Conclusions. The distance of the early-type stars in Westerlund 2 is established to be in excellent agreement with the distance of 
WR 20a, indicating that WR 20a actually belongs to the cluster. Our best estimate of the cluster distance thus amounts to 8.0 + 1 .4 kpc. 
Despite the earlier spectral types, the currently known population of early-type stars in Westerlund 2 does not provide enough ionizing 
photons to account for the radio emission of the RCW49 complex. This suggests that there might still exist a number of embedded 
early O-stars in RCW49. 
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1. Introduction 



The open cluster Westerlund 2 (Westerlund 1961) lies in a 
blowout region of the giant H n complex RCW 49 (Rodgers et 
al.[l960). The stellar winds and ionizing radiation of the early- 
type stars in Westerlund 2 have evacuated the dust in the central 
part of RCW 49 and filled the cavity with very hot low density 
gas. Unfortunately, so far, our knowledge of the stellar content of 
Westerlund 2 was rather fragmentary. Photometric studies were 
performed by Moffat & Vogt ( 119751 1, Moffat et al. ( 119911 1 and 
Carraro & Munari (2004). To date, the only spectroscopic study 
of the cluster was presented by Moffat et al. ( 1991 ) who obtained 
low resolution spectra of six O-stars in Westerlund 2 and classi- 
fied them as 06-7: V objects (although, as indicated by the colon, 
these spectral classifications were rather uncertain). 

Recently, interest in the stellar population of Westerlund 2 
was triggered by two independent observational studies. On 
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the one hand, RCW 49 was observed with the Infrared Array 
Camera (IRAC) aboard Spitzer in the framework of the 
GLIMPSE legacy su rvey ( Whitney et al. 120041 Churchwell et 
al. 120041 Churchwell 2005). These observations revealed strong 
evidence for ongoing star formation activity and underlined the 
need for a detailed study of the properties of the early-type stars 
in the cluster core and their impact on the surrounding nebula. 
On the other hand, WR 20a - one of the two Wolf-Rayet stars 
in RCW 49 - was found to be a very massive eclipsing binary 
consisting of two WN6ha stars with individual masses of about 
80 M (Rauw et al. |2004l|20031 Bonanos et al. [2004b . Since pre- 
vious spectroscopic studies revealed only one 06: and six 07: 
stars, in the cluster core (Moffat et al. 1991), the existence of a 
pair of stars that massive was somewhat unexpected. These re- 
sults therefore called for a re-investigation of the population of 
early-type stars in Westerlund 2. 



2. Observations 

2.1. Photometry 

Photometric monitoring of the core of Westerlund 2 was per- 
formed in service mode with the ANDICAM instrument at the 
1.3 m telescope at CTIO operated by the SMARTSQ consor- 
tium. The primary goal of this campaign was the study of the 
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Fig. 1. Left panel: average image (combining the best B and V filter observations, i.e. those with a seeing of less than 1.25 arcsec) 
of the region around Westerlund 2 as observed with the ANDICAM instrument. Right panel: zoom on the central part of the cluster. 
The labels correspond to the numbering scheme introduced by Moffat et al. ( 1991 ) and identify those stars for which we obtained 
spectra. Note that MSP 158 and 218 are late-type foreground stars that are unrelated to the cluster (see Sect. 13. II ). The boxes indicate 
the three newly discovered eclipsing binary systems discussed below. 



light curves of early-type stars and more specifically of the very 
massive eclipsing binary WR 20a. Whenever possible, the clus- 
ter was observed twice per night (separated by several hours) 
through the B and V filters of the ANDICAM instrument. Each 
pointing consisted of integration times of 2 x 200 and 2 X 75 s 
for the B and V filters respectively. The detector was a Fairchild 
447 CCD of 2048 x 2048 pixels. In the standard mode of the 
instrument the pixels are binned by 2 x 2. Each binned pixel cor- 
responds to 0.37 "on the sky and the total field of view of the 
instrument is about 6x6 arcmin 2 , well suited for monitoring si- 
multaneously most of the early-type stars in the central part of 
Westerlund 2 (see Fig.Q]). 

Our campaign covered 45 consecutive nights from 3 
December 2004 to 17 January 2005. The median seeing of our 
data is 1.4 arcsec. During photometric nights, observations of 
Landolt ( 1992) standard star fields were also performed. All the 
data were pipeline processed (bias subtraction and division by a 
normalized flat field) by the SMARTS consortium at Yale before 
they were delivered to us. 

The photometric reduction was done in three steps. First, in- 
strumental magnitudes were extracted with DAOPHOT (Stetson 
119871 ) under irai0 using aperture and PSF-fitting techniques. To 
avoid possible flat-fielding errors (Manfroid |1995l Manfroid et 
al. [2001 ) the atmospheric extinction and zero-point were deter- 
mined for the best photometric nights, using only the central part 
of the fields, by a multi-night, multi-star, multi-filter method as 
described by Manfroid (1993). The implementation of this re- 
duction procedure allows the construction of a consistent natural 
system, which contains the extra-atmospheric instrumental mag- 
nitudes of all constant stars included in the computation proce- 
dure. Stars from all observed fields were included in the build- 
ing of the natural system. To prevent adverse PSF effects, this 
procedure was carried out with the widest aperture photome- 
try and restricted to relatively bright, isolated stars. The exter- 
nal calibration is based on two Landolt (1992) standard fields, 
TPhe and Rubin 149, revisited by Stetson 0. These fields were 
observed at air masses ranging between 1.0 and 2.7. The rms 
error on the extinction coefficients of each night is 0.018. The 
error on the zero points is 0.028 mag, whilst the rms deviation 
of the mean magnitudes of the standard stars is slightly less than 
0.003 mag. Finally, an extended set of constant stars was built to 
provide a calibration for all the frames, including those of rel- 
atively poor quality, and over the whole field. This procedure 
(Manfroid 1995) yields robust differential data. As the centering 
of the fields was quite consistent over the whole data set, inaccu- 
racies in the overall flat-field calibration could not be evaluated, 
but they can be tolerated since they only consist of a systematic 
space-dependent shift and do not affect the variable component 
of the light curves. 

Comparing our PSF photometry to the values reported by 
Moffat et al. £1991) - calibrated by the photomultiplier aperture 

2 iraf is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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photometry of Moffat & Vogt ( 119751 ) - for 59 objects in common 
(excluding a few outliers) and Carraro & Munari d2004| ) for 82 
objects in common, we obtain: 

(B - V) MS p = (0.860 + 0.024)(fi - V) + (0.089 + 0.037) 
(B - V) CM = (0.958 + 0.016)(B - V) + (0.025 + 0.020) 
V MSP - V = (0.052 ± 0.029)(B - V) - (0.086 ± 0.044) 
V CM - V = (0.062 + 0.024)(B- V) - (0.183 + 0.029) 
where magnitudes and colours without an index refer to our 
data. As noted by Carraro & Munari (2004), the colour effects 
in the Moffat et al. (119911 ) data are large. Whilst our V zero- 
point agrees with the one of Moffat et al., the discrepancy with 
the zero-point of Carraro & Munari cannot be easily understood 
since most of the 82 stars in common are not affected by blends. 

2.2. Spectroscopy 

For the spectroscopic analysis, we adopt the MSP numbering 
convention introduced by Moffat et al. ( 119911 see also Fig.[T]). 
Medium-resolution long-slit spectroscopy of fourteen of the 
brightest stars in the central part of the Westerlund 2 cluster was 
obtained with the EMMI instrument at the ESO New Technology 
Telescope (NTT) on 27 March 2005. We excluded star MSP 91 
for which Moffat et al. dl991| ) obtained a spectral type GO V-III 
as well as MSP 273, another suspected non-member which lies 
rather far away from the cluster core. All seven O-type stars ob- 
served by Moffat et al. ( 1991 ) were included among our targets. 
Finally, MSP 168 was included in our target list, but could not 
be observed due to a lack of time. 

The red arm of EMMI was used in the REMD mode with 
grating #6 (1200 grooves mtrr 1 ). The useful wavelength domain 
ranged from 4200 to 4950 A. The spectral resolution of the in- 
strument configuration as measured from the FWHM of the lines in 
the ThAr calibration exposures was about 1 A. Exposure times 
ranged from 20 min to 1 hour. The spectra were bias-subtracted, 
flat-fielded, sky-subtracted, wavelength-calibrated and normal- 
ized using the midas software developed at ESO. The S/N ra- 
tios of the calibrated spectra ranged between 60 and 150 for the 
faintest and brightest sources respectively, with an average of 
~ 100 for a V = 14.5 star. Since the observations were taken 
only two days after full moon, the level of the background sky 
was rather high. Although we extracted the sky spectrum as close 
as possible to the source, some residual sky features are seen in 
the final spectra of the faintest source that we observed under 
somewhat poorer conditions (MSP 263). 

3. Spectral classifications 

3.1. Late-type foreground stars 

Two stars (MSP 158 and 218) have spectral types significantly 
later than B0 (see Fig. [2]). In order to establish the MK spectral 
types of these two objects, we make use of the digital spectral 
classification atlas compiled by R.O. Gray and available on the 
wetQ 

4 http://nedwww.ipac.caltech.edu/level5/Gray/frames.html 
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The spectrum of MSP218 (10:24:07.7, -57:45:41, J2000) is 
dominated by strong and broad H i Balmer lines. Weak absorp- 
tion lines of He i ^4388, 4471, Mgn A4481 and Sim A4552 
are also seen. The Hei /14471/Mgn ^1 448 1 intensity ratio in- 
dicates a spectral type B8-A0. The appearance of some weak 
metallic lines favours an AO classification, whilst the width of 
the Balmer lines suggests a giant luminosity class. We there- 
fore adopt an AO III spectral type. Our photometric data yield 
V = 14.13 + 0.03 and B - V = 0.34 + 0.04. If the intrinsic 
colours and absolute magnitudes of MSP 218 are those of a typ- 
ical AO III star (Schmidt-Kaler [19821 , then Ay ~ 1.15 and the 
star must be located at a distance of about 4.0 kpc. 

The spectrum of MSP 158 (10:23:59.5, -57:45:23, J2000) 
shows a prominent G-band as well as a strong Ca i A 4227 line in- 
dicating a late (G-K) spectral type. The intensity ratios between 
the Fei A 4325 and Hy lines as well as the ratio between Cri 
A 4254 and the neighbouring Fei /LI 4250, 4260 lines suggest a 
spectral type G8. The Yn /14376/Fei ,14383 line ratio and the 
CN bandhead at 4216 A indicate a giant luminosity class. We 
thus assign a G8III spectral type. From our photometric data, 
we infer V = 13.03 + 0.03 and B -V - 1.09 + 0.05. Assuming 
magnitudes and colours typical for a G8 III star (Schmidt-Kaler 
1982 ), we derive a reddening of Ay - 0.47 and a distance of 
about 2.25 kpc. 

Both stars are foreground objects unrelated to the 
Westerlund 2 cluster, as suspected by Moffat et al. d 1 99 It based 
on their UBV photometry. This is also confirmed by the absence 
of the strong diffuse interstellar band at 4430 A which is promi- 
nent in the spectra of the O-type members of the cluster (see 
Fig. [3}. As we will show below, the O-stars of Westerlund 2 are 
subject to substantially larger reddening (Ay from about 4.8 to 
5.8 mag) and their distance is significantly larger. Comparing the 
reddening of the G8 III, AO III and the O-type stars (see below), 
we note that the Av of the cluster members exceeds the value 
estimated from a crude empirical linear Ay(d) relation derived 
from the foreground stars. Though such a linear relation pro- 
vides only a first approximation, the result suggests that much of 
the reddening that affects the cluster could be local (e.g. due to 
material from the molecular cloud out of which the cluster has 
formed). 



3.2. Early-type stars 

Twelve of our targets display spectra typical of O-type stars 
(see Fig. 0). The spectral types were derived from the ratio 
of the equivalent widths of the Hei ,14471 and Hen A 4542 
classification lines by means of the criteria first introduced 
by Conti J1973I see also Mathys |1988l l. The results are dis- 
played in TableQ] For star MSP 203, we were only able to de- 
rive an upper limit on the equivalent width of the He i A 447 1 
line. All stars have rather strong Hen A 4686 absorption lines 
(equivalent widths in the range 0.4 to 1 . 1 A, the extremes cor- 
responding to MSP 167 and 263 respectively) suggesting they 
are main-sequence or giant objects. Mathys (1988) proposed 
EW = 0.56 A as the demarcation between weak and strong He n 
A 4686 absorption corresponding to luminosity classes III and V 
respectively. Using this criterion, we have assigned the luminos- 
ity classes in TableQ] 

We then compared our spectra to the digital atlas of Walborn 
& Fitzpatrick dT990l ) and Walborn et al. ( 120021 . The results are 
also indicated in TableQ] The classifications according to the two 
criteria agree to better than one subtype. We note that because 
the Hei A 4471 classification line is quite weak or broad in the 
spectra of MSP 167, 175 and 203, these objects could also be 
classified as 05 rather than 06 stars. Henceforth, we adopt the 
conservative hypothesis that these objects are actually 06 stars, 
but we will nevertheless consider the implications of a somewhat 
earlier classification for these stars when discussing the ioniza- 
tion of RCW 49. 

Three stars (MSP 18, 182 and 183) display a strong Hen 
A 4686 absorption accompanied by weak N m AA 4634-40 emis- 
sion, which leads to an ((f)) classification. The earliest star in 
our sample is MSP 183 which displays prominent N v /LI 4604, 
4620 absorption lines. Except for the weaker N m AA 4634-4640 
lines, the spectrum of MSP 183 very much resembles that of 
HD 64568 (03 V((f )), Walborn et al. [20Q2l We therefore clas- 
sify MSP 183 as 03V((f)). 



Fig. 3. Medium resolution spectra of O stars in Westerlund 2. 
Note the very strong diffuse interstellar band (DIB) at 4430A. 
Some residual nebular HB line emission affects the spectra of 
MSP 151, 182 and 199. 
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Fig. 2. Medium resolution spectra of the stars MSP 158 and 218 
in the field of Westerlund 2. Both stars are likely foreground ob- 
jects unrelated to the cluster. 



On average the spectral types that we determine are signifi- 
cantly earlier than those inferred by Moffat et al. (I1991I ). The im- 
plications of these revised classifications are discussed in Sect.[5] 
Uzpen et al. ([2005 ) recently inferred an 04 V((f)) spectral type 
for MSP 18. Our data have a superior resolution (1 A) than their 
observations (4.2 A) and are thus better suited to measure the 
strength of the weak He i classification lines. We therefore favor 
the 05 V((f)) spectral classification for this star. 

Though our data set is obviously not ideal for this purpose, 
we note that our spectra do not reveal strong evidence for mul- 
tiplicity. None of the stars has an obvious composite spectrum 
nor displays assymmetric, very broad or double lines that could 
make it a spectroscopic binary candidate. A higher resolution, 
high S/N spectroscopic monitoring of these stars is however 
needed to confirm or infirm this impression. 

In TableQ] we list the observed V magnitudes and the B - V 
colours for the spectroscopically identified O-type stars in our 
sample. For stars that are well isolated on the sky, the aperture 
and PSF photometry agree very well with each other and we 
adopted the magnitude obtained with an aperture of 8 pixels (i.e. 
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Table 1. Spectral classifications of the O-type stars in Westerlund 2. The first and fourth columns yield the number of the star in 
the numbering scheme of Moffat et al. ( 119911 ) as well as the spectral type proposed by these authors. The fifth and sixth columns 
provide our newly derived spectral types, determined using the Conti ( 1973) + Mathys ( 1988) and Walborn & Fitzpatrick ( 1990) 
classification criteria. Columns 7, 8 and 9 respectively yield our mean V magnitude, the B—V colour index and the extinction derived 
assuming the intrinsic colours from Martins & Plez (2006 ) and Ry — 3.1. The individual spectro-photometric distance moduli based 
on the calibration of Martins & Plez (2006 ) are listed in column 10. Column 11 provides the effective temperature as derived from 
our spectral classification and the calibration of Martins et al. (2005a). Finally, the bolometric luminosities given in the last column 
were computed for an adopted distance modulus of DM = 14.52 ± 0.38 (d = (8.0 + 1.4)kpc, see Sect. 15. 1 b . 



Star* 


RA 


DEC 




Spectral type 


V 


B-V 


A v 


DM 


T c ff 


(K) 


logL bo i/L G 




(J2000) 


(J2000) 


MSP 


This work 
























Conti crit. 


W&F crit. 
















MSP 18 


10:24:02.4 


-57:44:35 


07:V 


05.5 V-III 


05v-m((f)) 


12.88 


1.23 


4.68 


13.74 


40100 


±800 


5.89 ±0.16 


MSP 151 


10:24:00.5 


-57:45:24 


06:V 


07III 


06III 


14.33 


1.36 


5.05 


14.97 


37100 


± 1000 


5.37 ±0.17 


MSP 157 


10:24:00.8 


-57:45:26 


06:V 


05.5 V 


06.5 V 


14.33 


1.31 


4.93 


14.39 


38900 


± 1000 


5.37 ±0.17 


MSP 167 


10:24:02.1 


-57:45:28 


07:V 


06III 


06III 


15.18 


1.42 


5.24 


15.71 


38200 


±500 


5.14 ±0.16 


MSP 171 


10:24:04.9 


-57:45:28 




05 V 


04-5 V 


14.47 


1.45 


5.36 


14.53 


41900 


± 1000 


5.58 ±0.17 


MSP 175 


10:24:01.4 


-57:45:30 




05.5 V-III 


06V-III 


13.98 


1.29 


4.87 


14.55 


39000 


±800 


5.49 ±0.16 


MSP 182 


10:23:56.2 


-57:45:30 




04V-III 


04v-ni((f)) 


14.45 


1.32 


4.96 


15.30 


42600 


± 1000 


5.45 ±0.16 


MSP 183 


10:24:02.4 


-57:45:31 


07:V 


04 V 


03V((f)) 


13.57 


1.50 


5.52 


13.77 


43900 


± 1000 


6.07 ±0.17 


MSP 188 


10:24:01.2 


-57:45:31 


07:V 


C4V-III 


04 V-III 


13.41 


1.37 


5.12 


14.10 


42600 


± 1000 


5.93 ±0.16 


MSP 199 


10:24:02.7 


-57:45:34 




03 V 


03-4 V 


14.37 


1.40 


5.21 


14.88 


43900 


± 1000 


5.61 ±0.17 


MSP 203 


10:24:02.3 


-57:45:35 


Q7:V 


06 V-III 


06 V-III 


13.31 


1.37 


5.10 


13.59 


38600 


±800 


5.84 ±0.16 


MSP 263 


10:24:01.5 


-57:45:57 




Q6.5V 


06 V 


15.01 


1.61 


5.84 


14.09 


38350 


±500 


5.45 ±0.16 



3.0") radius which has the smallest errors. On the other hand, 
for stars in the crowded core of the cluster, it was mandatory to 
use the PSF photometry instead. Comparing the observed colour 
indices with the intrinsic colours of Martins & Plez (2006), we 
have inferred the extinction Ay, adopting R v = 3.1. Assuming 
that these objects have absolute magnitudes typical of main- 
sequence or giant O-type stars (Martins & Plez 2006), we esti- 
mate an average distance modulus of 14.47 + 0.65 (1-cr standard 
deviation). If we exclude the three stars with the lowest (MSP 18, 
183 and 203) and the star with the largest (MSP 167) value of 
the distance modulus, we get on average DM = 14.60 + 0.43 
in excellent agreement with the value (14.5 + 0.3) of Moffat 
et al. (119911 ) based on their UBV photometry and ZAMS fit- 
ting. This suggests therefore a distance of (8.3 + 1.6) kpc to- 
wards Westerlund 2 in very good agreement with the distance 
estimate of ~ 8.0 kpc towards WR20a (see Rauw et al. 2005 
and Sect. 14. It . We note that adopting the intrinsic colours from 
Schmidt-Kaler ( 1982 ) and the absolute magnitudes from Martins 
et al. (12005 al ) would reduce the distance modulus by 0.23 mag 
(i.e. a reduction of the distance by 800 pc). 



4. Photometric variability 

4.1. WR20a 

As pointed out in the introduction, one of the initial motivations 
of our observing campaign was the study of the photometric 
eclipses of the extremely massive binary WR20a. Bonanos et 
al. (|2004) reported / band observations showing that the system 
displays indeed eclipses on the orbital period of 3.686 days. The 
Fourier periodogram (computed using the method of Heck et al. 
1985 , revised by Gosset et al. 2001) of our photometric time se- 
ries of this system displays the highest peak at v' { = 0.5426 d~' 
for both filters (B and V) and irrespective of whether we con- 
sider the aperture or PSF photometry (see e.g. Fig.|4]). The nat- 
ural width of the peaks in the Fourier periodogram of our time 
series is equal to 0.0222 d -1 . Assuming that the position of the 
highest peak can be determined with an accuracy of 5% of this 
width, we infer an uncertainty of 0.001 1 d" 1 on v',. 




12 3 





III 















12 3 



F'requency v (d ) 

Fig. 4. Periodogram showing the half amplitude of the best-fit 
sine wave (top panel) and spectral window (bottom panel) de- 
rived from the Fourier analysis of the photometric time series 
of WR 20a. The results are shown for the measurements in the 
B filter (8 pixels radius aperture photometry). For the V filter, 
the results are identical. The highest peak in the periodogram is 
located at 0.5426 d _1 . 



Since the two components of WR 20a have roughly equal 
temperature and brightness, the eclipses are of essentially equal 
depth and shape. Thus v\ actually corresponds to twice the or- 
bital frequency. Therefore, our data yield an orbital period of 
(3.6859 ± 0.0075) days, in very good agreement with the value 
(3.686 + 0.010 days) found by Bonanos et al. ( 120041 . The whole 
periodogram (Fig.|4| can be interpreted on the basis of this sole 
periodicity. Indeed, no peak is visible at v ~ 0.27 d -1 because 
most of the power is contained in the v\ harmonic. The higher 
harmonic at v ~ 1.08 d -1 is induced by the shape of the eclipses 
whereas the harmonics at v ~ 0.82 d -1 (not to be confused with 
v = 0.93 d -1 , the one-day alias of 1.08 d -1 ) is related to the slight 
difference in eclipse depths. 
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Folding the light curve of WR 20a with the derived period, 
we have evaluated the magnitudes of the system at phases out- 
side of the eclipses: B = 14.999 ± 0.021, V = 13.416 ± 0.024 
and / = 10.658 + 0.015 (the latter value was computed from the 
Bonanos et al. data). Our estimate of the period can be further 
refined by comparing our light curves with that of Bonanos et 
al. Subtracting the above outside-eclipse magnitudes from the 
observed light curves, we can combine all the data and per- 
form a Fourier analysis. In this way, we find that the highest 
peak in the periodogram has v x = (0.54277 + 0.00019) d" 1 . 
Therefore, the newly determined value of the orbital period is 
(3.68475 + 0.00129) days. The value of the uncertainty is prob- 
ably too conservative: comparing the time of primary minimum 
determined by Bonanos et al. (120041 1 with our best estimate of 
the primary minimum results in a more realistic estimate of the 
uncertainty on the orbital period of 0.00020 days. 

The light curves of WR 20a folded with the above period 
and adopting a time of primary minimum of HJD 2453330.930 
are displayed in Fig. [5] We have fitted our B and V light curve 
along with the / band light curve from Bonanos et al. ((2004) 
using the nightfali0 software developed and maintained by R. 
Wichmann, M. Kuster & R Risse. We fixed the mass ratio at 
the value (m sec /m pr i m = 0.99) derived from the radial velocity 
curve. There exist a number of combinations of the model pa- 
rameters that fit the light curve but do not provide a physically 
self-consistent description of the system. For instance, if the fill- 
ing factors (defined as the ratio of the polar radius of the star 
to the polar radius of the Roche lobe) of both components are 
allowed to differ in the fit of the light curve, the luminosities of 
the primary and secondary star are found to differ by more than 
50% which is clearly at odds with the fact that the two stars have 
identical spectral types and spectral features of identical strength 
(Rauw et al. [2004 2005). We have therefore constrained some 
of the fitting parameters. For instance, we searched for models 
which have primary and secondary temperatures in agreement 
with the value (43000 + 2000 K) derived from an analysis with 
a model atmosphere code (see Rauw et al. 120051 1. Following 
the same strategy as Bonanos et al. (2004), we have assumed 
that both stars have identical (Roche lobe) filling factors. For a 
grid of values of the filling factors and the primary (resp. sec- 
ondary) temperature, we iterated the nightfall code for the best 
fit inclination and secondary (resp. primary) temperature. The 
best constrained model parameter is the orbital inclination which 
quickly converges to i = (74.5 + 1.0)° in excellent agreement 
with the value derived by Bonanos et al. (|2004). Moreover, the 
shape of the wings of the eclipses and the light variations out- 
side of the eclipses require that the filling factor be between 0.89 
and 0.92 with a preferred value of 0.91. The absolute tempera- 
tures of the components cannot be constrained by the analysis of 
the light curve. In principle, the depths of the eclipses constrain 
the relative surface brightnesses and hence the relative tempera- 
tures. Under the above assumptions of equal radii, we find that 
the primary should be hotter than the secondary by about 2000 
- 3000 K (assuming a primary temperature in the range found 
from the spectroscopic analysis), although models with identi- 
cal temperatures for both components are of almost equivalent 
quality. The best fit light curves for a filling factor of 0.91 for 
both stars, an orbital inclination of 74.5°, and temperatures of 
T efi .p = 43000 K and T effjS = 40500 K are illustrated in Fig.0 



5 http://www.hs.uni-hamburg.de/DE/Tns/per/Wichmann/Nightfall.html 
Note that the nightfall code does not include stellar wind effects, but 
given the proximity of the stars in WR20a, the light curve should be 
dominated by the geometrical effects. 



Fig. 5. Light curve of WR20a in th e B, V and / (the lat- 
ter data from Bonanos et al. 2004) filters folded with the 
3.68475 day period and adopting a time of primary minimum of 
HJD 2453330.930. The B and V data correspond to the photom- 
etry obtained through PSF fitting and are given in Table 2 (avail- 
able in electronic form at the CDS). Filled symbols stand for data 
with a FWHM of the PSF of less than 4.0 pixels, whilst the open 
symbols yield data taken under poorer seeing conditions. The 
continuous curve represents the best fit theoretical light curve 
determined with the nightfall code for a filling factor of 0.91 
for both stars, an orbital inclination of 74.5°, and temperatures 
of T eff , p = 43000 K and T eff , s = 40500 K. 



The parameters of the above solution correspond to an av- 
erage stellar radius of (18.7 + 0.9) R , roughly 5% smaller 
than the values found by Bonanos et al. (2004). If the pri- 
mary temperature amounts to (43000 + 2000) K as found by 
Rauw et al. (2005), then its bolometric luminosity would be 
log L/Lq = 6.03 + 0.09, whilst for the secondary we obtain 
log L/L G = 5.93+0.10. For the entire binary system, we thus find 
log L/L Q = 6.28 + 0.10. With the bolometric correction of -3.91 
(Rauw et al. [2005), we therefore infer an absolute V magnitude 
of My = -7.04+0.25. The out of eclipse magnitudes listed above 
yield B - V - 1.58 ± 0.03, hence E(B - V) = 1.91 ± 0.03 (eval- 
uated using the intrinsic colour index from Rauw et al. 2005) 
and A v = 5.93 + 0.09. This leads to DM = 14.52 + 0.27 cor- 
responding to d = (8.0 ± 1.0) kpc. This result concurs with the 
estimate of Rauw et al. (2005, d — (7.9 + 0.6) kpc), and agrees 
perfectly well with the estimate of the distance of the O-type star 
population presented above. We will come back to this issue in 
Sect.0 

Although the periodogram of Fig. [4] can be fully explained 
by the orbital period and its harmonics, we decided to search 
for additional variations, in order to extract the maximum in- 
formation from the data. Moreover, it is reasonable to think 
that such massive stars might display additional variations such 
as pulsations (see Sterken & Breysacher 1997, and Fig. 3 of 
Pamyatnykh 1999). Consequently, we have whitened the B, V 
photometry time series for the binary variations by removing the 
orbital period and its harmonics. The time series of the residuals 
have again been Fourier analysed. No peak exceeding a semi- 
amplitude of 0.008 mag has been found, indicating that no other 
periodicity is present in the photometric variations of WR 20a at 
this level. 

4.2. O-type stars 

Using our photometric time series, we have searched for vari- 
ability among the objects identified as early-type O-stars in our 
spectroscopic campaign. For targets that are rather well isolated, 
the aperture photometry yields 1-cr dispersions on the B and V 
magnitudes of order 0.005 to 0.010 mag. A Fourier analysis of 
our time series, excluding the data points with the poorest see- 
ing reveals no significant peak in the periodogram. The highest 
peak is usually found near v - 1 d~' which is actually due to 
the interplay between the sampling of our time series and a few 
deviating data points. 

For objects in the crowded region of the cluster core, we find 
that the variations of the aperture photometry time series reflect 
the variations of the seeing during the observing campaign as 
could be expected. Excluding those points with the poorest ob- 
serving conditions from the aperture photometry time series, we 
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Fig. 6. L<?/f: light curves of MSP 223 as obtained from our B and V PSF photometry time series folded with a period of 29.21 1 days 
and Tq = 2453369.85. The filled and open symbols stand for data obtained with a seeing of less than and more than 4.0 pixels 
(fwhm) respectively. Middle: same as left panel, but for MSP 96, a period of 1.0728 days and To - 2453370.6973. Right: same as 
left panel, but for MSP 44, a period of 5. 176 days and T = 2453371.75. 



find that there is no photometric variability exceeding the l-cr 
level of 0.010 mag. The Fourier analysis yields the same conclu- 
sions as for the isolated objects. Although the PSF photometry 
provides more accurate absolute photometry, PSF photometry 
time series usually yield somewhat larger 1-cr dispersions. We 
conclude that none of the early O-type stars in the cluster core 
displays a significant variability and in particular, we find no ev- 
idence for photometric eclipses. 



4.3. Other variables in the field of view 

We have searched our photometric time series for additional 
variables. The light curves of those objects that display a pho- 
tometric variability exceeding significantly the level of photon 
noise in the B and V band were extracted and were analysed 
with the Fourier technique over the frequency range between 0.0 
and 5.0 d -1 . Only data points obtained when the seeing was bet- 
ter than 4 pixels (i.e. 1.5") were included in this process. For 
most objects, the Fourier analysis yields the highest peaks at fre- 
quencies near or an integer number of day -1 . These features 
stem from the combination of the actual sampling of the time 
series with a few deviating data points and do not reflect a gen- 
uine periodic behaviour. The only exceptions are the three stars 
discussed below. 

MSP 223 (10:23:59.2 -57:45:40, J2000) is an apparently red 
star (V = 15.75, B -V - 1.50) that shows drops in its light 
curve that are consistently seen in the B and V filter and occur 
on a recurrence time of 14.61 days. The variations are reminis- 
cent of an eclipsing binary. In this case, the system either lacks 
a secondary eclipse (maybe as a result of an eccentric orbit or 
of an insufficient sampling of the light curve near that phase) or 
the true orbital period could be 29.21 days. We folded the light 
curve with this latter period and adopting a time of primary min- 
imum of HJD 2453369.85. The eclipses are rather narrow and 
the secondary eclipse could be slightly shallower than the pri- 
mary one. The system presents some slight slow variations out- 
side the eclipses (see Fig.|6]l. Though the photometric data are 
not sufficient to provide a full solution, test calculations with the 
nightfall code suggest that the system consists of a pair of stars 
of comparable surface brightness that are well inside their Roche 
lobe and are seen under a high inclination angle. The variations 
outside eclipse cannot be explained by ellipsoidal variations and 
rather suggest the existence of a non-uniform surface brightness 



for at least one component. MSP 223 could thus be a kind of 
magnetically active binary. 

The periodogram of the photometry of MSP 96 (10:24:01.0 
-57:45:06, J2000; V = 16.17, B - V = 1.27) displays a promi- 
nent peak at 1.864day _1 , although its one-day aliases are only 
slightly weaker. Again the light curve suggests an eclipsing bi- 
nary and the true period is probably 1 .0728 day. The eclipses are 
rather broad (see Fig.|6j and the ellipsoidal variations suggest 
that at least one of the two stars must be close to filling its Roche 
lobe. Also, the different depths of the eclipses indicate that one 
of the stars must be significantly hotter than its companion. 

The periodograms of MSP 44 (10:24:00.5 -57:44:45, J2000; 
V = 15.67, B - V = 1.25) display their highest peak at 
0.1932day~', although the 1 - v alias is almost equally strong. 
The light curve folded with a period of 5. 176 days and To = 
2453371.75 displays a strong minimum at phase 0.0 as well as a 
much weaker secondary minimum at phase 0.5. Outside these 
rather narrow minima the light curve is not flat, but displays 
a strong modulation (see Fig. [6}. As for MSP 223, the narrow 
eclipses are the signature of a pair of relatively compact stars 
seen under high inclination, whilst the strong modulation of the 
flux outside eclipse requires a non-uniform surface brightness 
(spots or reflection effects) probably for both components. The 
light curve of MSP 44 is reminiscent of those of RS C Vn sys- 
tems therefore making this a candidate for a magnetically ac- 
tive binary system. Alternatively, MSP 44 could be a detached 
binary with components that differ quite substantially in temper- 
ature. The seemingly symmetric variability outside of the eclipse 
would then result from a strong reflection effect or surface heat- 
ing. 

To further constrain the nature of MSP 44, 96 and 223, we 
examine the (B - V, V) colour-magnitude diagram of the star^]in 
the ANDICAM field of view around Westerlund 2 (see Fig.|7J. 
Figure|7] clearly reveals two populations of stars: the foreground 
field stars which have relatively blue colours and the cluster stars 
which have an E(B - V) colour excess between 1.51 and 1.88. 
The O-type stars and WR20a fall in the region of the cluster 
main-sequence (see also Sect. |5.2| i. Two of the three foreground 
stars MSP 91 and 218 clearly occupy the locus of field stars. 
The third one (MSP 158, G8III) falls close to the cluster main 
sequence, suggesting that the locus of the cluster stars might 
be polluted by cluster non-members. The three eclipsing bina- 

6 Tables 3 and 4 with the light curves of the eclipsing binaries and 
the mean (B - V, V) data are provided at the CDS. 
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Fig. 7. Colour-magnitude diagram of the stars in the ANDICAM 
field of view around Westerlund 2. The filled circles and the open 
triangle correspond to the twelve O-type stars observed in spec- 
troscopy (see Table Q]) and the WR 20a binary, respectively. The 
dashed lines correspond to the main-sequence between spec- 
tral types 03 and B2 for a distance modulus of 14.52 and red- 
dened with Ay = 4.68 (leftmost curve) or 5.84 (rightmost curve) 
and assuming Ry = 3.1. The open circles yield the late-type 
foreground objects MSP 91, 158 and 218 whilst the filled trian- 
gles indicate the three eclipsing binaries discussed in Sect. 14. 31 
Finally, the two squares correspond to the two non-periodic vari- 
ables. 



ries fall on the locus of the cluster main-sequence. Moffat et al. 
d 1 99 It actually labelled all three stars as probable OB members 
of Westerlund 2 based on their UBV photometry. In this case, 
the primary components of these binaries would be late O or 
early B-type stars which would contradict our suggestion of two 
of these systems being magnetically active late-type binaries. A 
spectroscopic follow up of these objects will not only help clar- 
ify their nature, but will also bring valuable constraints on the 
general cluster properties. 

Finally, Kolmogorov-Smirnov and x 1 tests were also per- 
formed to identify non-periodic variables. Two more stars 
emerged from these tests and both were found to display a 
strong apparently non-periodic variability. These objects are lo- 
cated at 10:23:58 -57:47.7 and 10:24:00 -57:44.9 (J2000) and 
have magnitudes and colours of V = 20.01, B — V = 1.83 and 
V = 18.93, B — V — 1.88, respectively. Both objects are about 
1 mag fainter (in B and V) at the beginning of our campaign 
than later and these 'minima' last about 14 and 18 days respec- 
tively. No other event of this type is seen over the rest of the 
campaign. These objects are either irregular or long-term (pe- 
riod > 45 days) variables or eclipsing binaries with periods very 
close to 1 or 2 days. 



5. The link between Westerlund 2, WR20a and 
RCW49 

5.1. The distance 

The distance towards Westerlund 2 and RCW 49 has been a very 
controversial issue (see e.g. the discussions in Churchwell et 
al. I2004 and Rauw et al. I2005l l and distance estimates for the 
cluster range from about 2.5 kpc to 8kpc. The line of sight to- 
wards RCW 49 is roughly tangential to the Carina arm, lead- 
ing to a wide range of possible distances. Using a galactic ro- 
tation curve along with H i line velocity measurements, Russeil 
(2003) recently determined a kinematic distance of 4.7tS'| kpc 
towards RCW 49. Russeil cautioned however that velocity devi- 
ations from the circular Galactic rotation model exist all along 
the Carina spiral arm. 

Based on spectro-photometric parallaxes of the early O-type 
stars, we have estimated d = (8.3 ± 1.6) kpc, whilst the anal- 
ysis of the light curve of the eclipsing binary WR20a yields 
d = (8.0 + 1.0) kpc. The weighted mean of these two distance 
estimates amounts to d = (8.0 + 1.4) kpc. These results clearly 
rule out a distance below ~ 6 kpc. Moreover, the excellent agree- 
ment between our two distance determinations indicates that 
WR 20a very likely belongs to Westerlund 2. Various effects that 
could bias our distance determinations are briefly discussed in 
AppendixlAl 

Based on Spitzer data, Uzpen et al. (j2005) found that 
MSP 18 displays an IR excess which they attributed to ther- 
mal bremsstrahlung from the stellar wind. These authors fitted 
the V, Spitzer and 2MASS photometry with a Kurucz ATLAS9 
model. In this procedure, the scaling factor of the surface flux 
depends on the angular diameter (i.e. on the ratio of the stellar 
radius and the distance) of the star. In this way, they inferred 
a spectro-photometric distance of 3233*f«pc and an extinction 
of Ay = 5.63*Qji. These authors suggest therefore that the ac- 
tual distance of Westerlund 2 is 3.2 kpc. The extinction is signif- 
icantly larger than our value (4.68), whilst their distance is less 
than 50% of our best estimate of the cluster distance. To under- 
stand the origin of this discrepancy, we fitted our B, V magni- 
tudes along with the U magnitude from Moffat et al. ( 1991) and 
the 2MASS photometry using a simple (static, plane-parallel at- 
mosphere in LTE) model atmosphere codfl Test fits including 
the 2MASS data suggest that the large Ay value of Uzpen et al. 
(2005) is likely due to the near-IR domain already being affected 
by the IR excess emission. Freezing the reddening at the value 
of E(B — V) = 1 .56 found above, the temperature at 40 kK com- 
patible with the optical spectrum and fitting only the UBV data 
clearly confirms this idea. The angular diameter derived from the 
latter fit amounts to 6 x 10~ n rad. Assuming a typical radius of 
12.0 R© for the 05 V-III star then yields a distance of 9.0 kpc. 
We thus conclude that the discrepancy between the Uzpen et al. 
(120051 1 result and our estimate stems from the near-IR photome- 
try which is already affected by excess emission. 

In summary, we find that the distance of WR20a and 
Westerlund 2 is probably (8.0 + 1.4) kpc. 



5.2. The Hertzsprung-Russell diagram of Westerlund 2 

To position the O-type stars in the Hertzsprung-Russell diagram, 
we have adopted the 'observed' spectral-type - effective tem- 
perature calibration and bolometric corrections of Martins et al. 



7 The software, developed by S. Jeffery at Armagh Observatory, is 
available at http://star.arm.ac.uk/~csj/software.store 



G. Rauw et al.: Early-type stars in the core of the young open cluster Westerlund 2 



8 

(2005a) and Martins & Plez (2006), respectively. The observed 
V magnitude and Ay extinction of each star were taken from 
Table Q] The bolometric luminosities are quoted in the last col- 
umn of Table Q] and were evaluated accounting for the uncertain- 
ties in spectral type and luminosity class and assuming a distance 
of (8.0+ 1.4)kpc. 

In Fig. [8] we compare the location of the O-stars in the 
Hertzsprung-Russell diagram with the evolutionary tracks and 
isochrones from Schaller et al. (1992]) for solar metallicity. 
Keeping in mind the possibility that some of the luminosities 
might be overestimated due to binarity, we note that the O-stars 
that we have observed should have masses ranging from about 
30 to 85 M and ages of less than -2.5 Myr. 





4.8 4.7 4.6 4.5 4.4 4.3 

l0 g T eff 

Fig. 8. Hertzsprung-Russell diagram of the O-type stars in 
Westerlund 2. Evolutionary tracks are from Schaller et al. 
( 1 19921 ). The dashed line yields the 2.5 Myr isochrone corre- 
sponding to an upper limit on the age of the O-stars. The filled 
triangle indicates the position of the WN6ha stars of WR 20a. 



5.3. The ionization of RCW49 

Several authors have presented radio observations of RCW49. 
The continuum radio emission of H n regions is free-free emis- 
sion that can, in principle, be used to evaluate the minimum num- 
ber of Lyman continuum photons. In the case of RCW 49, the ob- 
served radio luminosity (evaluated for a distance near 8 kpc) was 
usually found to be too large compared to the Lyman ionizing 
flux from the known O-star population. For instance, Whiteoak 
& Uchida ( [1997] ) evaluated a flux density of 210 Jy at 0.843 GHz 
for the entire RCW 49 complex, whilst the flux density of the 
shell around Westerlund 2 was found equal to 110 Jy. Based 
on these numbers and the census of O-type stars of Moffat et 
al. (119911 ), Whiteoak & Uchida (119971 ) further argued that the 
Lyman continuum flux of the early-type stars could only explain 
the observed radio emission of RCW 49 if the distance were of 
order 2.3 kpc. 

In this context, we emphasize that adopting a lower distance 
for the nebula and the cluster does not solve the issue of the 



missing ionizing photons. In fact, both, the radio flux and hence 
the required number of ionizing photons on the one hand and the 
bolometric luminosity of the stars on the other hand scale with 
the distance squared. Assuming a lower distance as was done by 
Whiteoak & Uchida d 19971 ) would thus imply a drastic under- 
luminosity of the stars compared to typical objects of the same 
spectral type. Therefore, the ionizing flux would also be reduced 
compared to the values for typical stars. 

For a distance of 8.0 kpc, the number of Lyman continuum 
photons required to account for the total 0.843 GHz flux of 
RCW 49 and the flux of the shell near Westerlund 2 as quoted 
by Whiteoak & Uchida ( [19971 ) would be 11.5 x 10 50 and 
6.0 x 10 50 s _1 , respectively. In a similar way, scaling the result 
of Churchwell et al. (120041) to our preferred distance estimate 
yields a minimum UV photon luminosity of 20.7 x 10 50 s _1 to 
account for the 5 GHz flux density of 335 Jy (as determined by 
Goss & Shaver |T9701 ). We note that lower determinations of the 
5 GHz flux can be found in the literature: 178.8 Jy (Smith et al. 
[1978] ) and 263.2 Jy (Kuchar & Clark [T997l ), which imply photon 
numbers in the range 11.1 - 16.3 x 10 50 s _1 . On average, we 
find (14.9 + 4.5) x 10 5() s _1 . These numbers actually represent a 
lower limit on the ionizing flux since they neglect the possible 
leakage of some UV radiation from the nebula (see Churchwell 
et al. 120041 ) and they were estimated for a dust- free nebula using 
the formalism of Mezger et al. dl974| ). Near- and mid-IR images 
of the RCW 49 complex reveal however the existence of dust 
in this nebula, although the stellar winds and radiations of the 
O-type stars have cleared the region immediately surrounding 
W esterlu nd 2 (e.g. Conti & Crowther 2004] a nd Ch urchwell et 
al. 120041 ). Using the formalism of Smith et al. ( 119781 ) to account 
for the impact of dust, we estimate that about 58% of the UV 
continuum flux is actually used to ionize the nebula. As a 
result, the total number of Lyman photons required to maintain 
the ionization of RCW 49 would have to be (25 .7 +7.9) x 10 50 s" 1 . 

With the newly determined spectral types and bolometric lu- 
minosities at hand, we have evaluated the number of ionizing 
photons produced by the O-type stars of Westerlund 2. To this 
aim, we have scaled the numbers of Lyman ionizing photons 
emitted per second (Qo) as tabulated by Martins et al. (f2005a) 
by the ratio of the above determined bolometric luminosity di- 
vided by the 'typical' theoretical bolometric luminosity. For a 
distance of 8.0 kpc, the total number of Lyman ionizing photons 
produced by the twelve O-type stars in our sample amounts to 
3.2 x 10 50 s _1 . Assuming that stars in the spectral range 06.5 - 
03 have masses between 28 and 60 M Q , and adopting a Salpeter 
IMF, we have estimated that Westerlund 2 should harbour about 
20 stars with spectral types in the range 09.5 - 07 (masses of 
15.5 - 28 M ). These late O-type stars would contribute about 
7 x 10 49 s~' (i.e. ~ 20% of the early-type O-stars) to the total 
number of Lyman ionizing photons. Adding the ionizing flux of 
the W R20a pair of WN6ha stars (1.45 x lO^s -1 , Rauw et al. 
120051 ) we arrive at a total of 5.4 x 10 50 s -1 , i.e. a fifth of the pho- 
ton flux required to explain the radio emission of RCW 49 (or 
slightly more than one third of the required flux if we neglect 
the impact of the dust). In Sect. 13. 21 we noted that three stars 
of our sample (MSP 167, 175 and 203) might actually be some- 
what earlier than their classification listed in TableQ] Assuming 
these objects to be of type 05, rather than 06, would increase 
the amount of ionizing photons by 1.5 x 10 49 s _1 , i.e. by about 
3%. This latter remark does certainly not alter our conclusion 
about the lack of a sufficient number of ionizing photons. 

Of course, our estimate provides only a lower limit on 
the actual production of ionizing photons inside the RCW 49 
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complex. For instance, we have not included WR20b in our 
account. This star was classified as WN7:h by Shara et al. 
dl9991 l. By analogy with the ionizing fluxes of WN6ha stars 
in NGC 3603 (Crowther & De ssart [T998i the components 
of WR20 a (Rauw et al. [2Q05l > and the models of Smith et 
al. (120021) . we estimate that log(Qo) < 50.0. Thus this star 
alone cannot solve the issue of the missing ionizing photons. 
Some of the early-type stars might still be hidden from our 
view by the material of their natal cloud. Indeed, observations 
indicate that star formation is ongoing in Westerlund 2. Whitney 
et al. (120041) reported Spitzer IRAC near-IR colours of the 
stars in RCW49 that are typical of a star formation region. 
They estimated the total number of young stellar objects to 
be ~ 7000. Star formation occurs preferentially within ~ 5' 
from the cluster core, where the stellar winds of the massive 
stars in Westerlund 2 have swept up gas and dust. This suggests 
that the ongoing star formation activity may be induced by 
the action of the early-type stars in Westerlund 2. The spectral 
energy distributions of two stars were fitted with models of 
B2 V and B5 V stars surrounded by accretion disks and a large 
envelope feeding the disk (Whitney et al. 120041 Churchwell 
120051 ). However, these stars are apparently not massive enough 
to contribute significantly to the ionization of the nebula. If the 
predictions of stellar atmosphere models are correct, then we 
conclude that we must still be missing an important fraction of 
the early-type stars in RCW 49. 

The influence of the early-type stars in the cluster core is 
not limited to the ionization of the nebula. The stellar winds of 
these stars also strongly impact on their surroundings. The sig- 
nature of this process has been detected as diffuse X-ray emis- 
sion with Chandra (Townsley et al. 120051 ) . This feature is proba- 
bly due to the thermalisation of the stellar winds either resulting 
from the collisions between the winds of individual early-type 
cluster members or from the shocks with the ambient interstel- 
lar medium. Townsley et al. (2005 ) found that this emission has 
temperatures between kT ~ 0.8 and 3.1keV. Assuming a dis- 
tance of 2.3 kpc, these authors inferred a total intrinsic (i.e. ab- 
sorption corrected) X-ray luminosity of 3 x 10 33 erg s _1 over the 
0.5 - 8.0 keV energy range. For a distance of 8.0 kpc, this X-ray 
luminosity becomes 3.6 x 10 34 ergs~'. From the temperatures 
and emission measures of the best-fit X-ray models proposed 
by Townsley et al. (2005), we estimate that a total thermal en- 
ergy of 6.5 - 20.4 x 10 49 erg (for a filling factor of 0.1 and 1.0 
respectively) is required to produce the observed diffuse X-ray 
emission. 

Considering a Salpeter IMF and the typical masses of O- 
stars as tabulated by Martins et al. (2005a), we estimate that 
Westerlund 2 should contain about 4500 M in the form of stars 
with masses between 1 and 120 M . We used the starburst99 
code (v5.0; Leitherer et al. 1 19991 Vazquez & Leitherer l2005l) to 
estimate the total mechanical energy injected through the stellar 
winds. For an age of 2Myr, no supernova explosion should have 
occured. However, the total wind luminosity should be about 
5.4 X 10 37 ergs _1 and the cumulative energy input from stellar 
winds would be about 3.6 x 10 51 erg. This exceeds the thermal 
energy that produces the X-ray emission by more than one or- 
der of magnitude. Since RCW 49 is a blister-like structure, the 
latter result is not surprising. Recent models accounting for the 
effects of clumping indicate that the mass loss rates of O-type 
stars might have been overestimated and are probably lower by 
a factor 3-4 (e.g. Martins et al. I2005bl Bouret et al. [2005 ). 
Lower mass loss rates would imply a reduction of the cumulative 
energy input from the stellar winds, although the latter number 



would still exceed the thermal energy responsible for the X-ray 
emission by at least a factor of a few. 

In summary, while the currently known population of O-type 
stars in Westerlund 2 does apparently not provide enough ioniz- 
ing photons to explain the radio emission of RCW 49, the dy- 
namical energy provided by the stellar winds of these stars seems 
sufficient to explain the observed diffuse X-ray emission. 

6. Summary and conclusions 

Our spectroscopic observations indicate that the core of 
Westerlund 2 harbours at least twelve O-type stars of spectral 
type earlier than 07, the earliest object being of type 03 V((f)). 
These spectral types are significantly earlier than previous classi- 
fications. The spectro-photometric distance of these O-type stars 
is found to be 8.3 + 1.6 kpc. This number is in excellent agree- 
ment with the result obtained from the study of the B and V 
light curve of the very massive Wolf-Rayet binary WR 20a and 
the distance based on ZAMS fitting of the UBV photometry 
of Moffat et al. d 1 99 It . The weighted distance of Westerlund 2 
and WR20a amounts to 8.0 ± 1.4kpc. This distance determi- 
nation implies that the radio emission of the giant Hn region 
RCW 49 must be larger than previously estimated. Despite the 
earlier spectral types inferred in this paper, the currently known 
population of early-type stars in Westerlund 2 provides only be- 
tween one fifth and one third of the ionizing photons required 
to account for the radio emission of RCW 49. Future observa- 
tions should aim at establishing a complete census of the O-type 
stars population in RCW 49, although our colour-magnitude dia- 
gram (Fig. [7]) indicates that at least within the field of view of the 
ANDICAM instrument, there are only a few candidates for early 
O-type stars other than those objects we have classified here. 

None of the early O-type stars in the cluster core displays a 
significant photometric variability on time scales from a day to 
a month and a half. None of these objects displays an obvious 
signature of multiplicity in our snapshot spectra. However, we 
have identified three eclipsing binaries (in addition to WR 20a) 
that are probably late O or early B-type stars belonging to the 
Westerlund 2 cluster. This suggests that there might in fact be 
many more (non-eclipsing) binaries. A spectroscopic follow-up 
of these systems as well as of some of the O-type stars that 
appear somewhat overluminous (MSP 18, 183 and 203) would 
bring precious information on their nature as well as on the gen- 
eral properties of the cluster. 

Finally, we emphasize that the results reported here indicate 
that WR 20a actually belongs to the Westerlund 2 cluster and fur- 
ther suggest that the cluster contains enough massive stars to ex- 
plain the location of the Wolf-Rayet binary through the effect 
of a dynamical interaction scenario as suggested by Rauw et al. 
(120051 

Appendix A: Uncertainties of the distance 
determinations 

The distance of WR 20a stems from the luminosity determined 
in the light curve analysis. Under the assumptions made in Sect. 
14.11 the radii of the stars in WR 20a are rather robust. The issue 
of the effective temperature is a bit trickier. This parameter is 
essentially determined by the spectroscopic analysis with a non- 
LTE model atmosphere code (Rauw et al. 2005). The strengths 
of the N in, N iv and N v features provide the most stringent con- 
straints. The model atmosphere code used by Rauw et al. (2005) 
assumes a spherically symmetric star. In a binary system such 
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as WR 20a, the stars are however not spherical. Their shape is 
rather set by the Roche potential. Moreover, the gravity varies 
over the stellar surface, being maximum near the poles and min- 
imum on the portion of the surface facing the companion star. 
As a result, gravity darkening renders the surface temperature 
non-uniform: the star is hotter near the poles and cooler near the 
equator. Adopting the fitting parameters of the light curve shown 
in Fig. [5] and assuming a gravity darkening exponent /3 — 1.0, 
we estimate that the temperature near the primary's poles should 
be about 1900K hotter than the average effective temperature 
(43000 K), while it would be lower by ~ 5300 K over the part 
of the stellar surface facing the secondary. Whilst these numbers 
may be overestimated, for instance because we neglected reflec- 
tion effects, they nevertheless suggest that part of the difficulties 
encountered when trying to fit simultaneously N m, N iv and N v 
lines might actually arise from a non-uniform surface temper- 
ature distribution in a close binary system. Concerning the dis- 
tance of WR 20a, we estimate that the effect of gravity darkening 
could produce a 4 - 5% uncertainty on the effective temperature 
which would then correspond to a 10% margin on the distance. 

Alternatively, our distance estimate of the early O-type stars 
in Westerlund 2 could be biased. In this respect, we note the 
rather large dispersion (0.65 mag) of our spectro-photometric 
distance moduli. The primary suspects to bias our result would 
be binary systems. However, our photometric monitoring cam- 
paign failed to reveal any eclipsing binary system among these 
O-type stars with a period between a few days and a few weeks. 
Furthermore, we do not find evidence for multiplicity in our 
spectroscopic data either, although these snapshot data are cer- 
tainly not suited to provide firm constraints on the number of 
spectroscopic binaries. Therefore, it seems unlikely that the ma- 
jority of the O-type stars that we have observed are binaries. If 
we nevertheless consider that some of them might be multiple 
objects, the smaller distance moduli of MSP 18, 183 and 203 
make them the best binary candidates. Discarding these three 
objects from the mean yields DM = 14.72 + 0.54 (1-cr standard 
deviation). Note that the dispersion about the mean remains quite 
large, though comparable to the cosmic scatter of absolute mag- 
nitudes based on the spectral types. 

Another possible reason for the dispersion of DM could be 
that some of the main-sequence O-stars are slightly evolved off 
the main sequence and are thus intrinsically brighter than 'typ- 
ical' main-sequence stars. Our luminosity class determination 
based on the strength of the He n A 4686 line suggests that some 
of the stars are indeed giants. Those stars that have an EW of 
the He n A 4686 line at the border-line between the two luminos- 
ity classes are labelled with a V-III luminosity class in TableQ] 
For these objects, we evaluated the distance by comparison with 
an interpolation of the main-sequence and giant absolute magni- 
tudes of Martins & Plez (2006). Some other objects could have 
intermediate evolutionary stages and might thus be responsible 
for part of the dispersion. 

Finally, as suggested in Sect. 13. II some part of the redden- 
ing that affects the stars in Westerlund 2 might be produced by 
material from the parental molecular cloud. Therefore, we can- 
not rule out a priori that, at least for this material, Ry differs 
from the canonical value of 3.1. In fact, in some clusters with 
a higher mean density interstellar medium, it has been argued 
that Ry could range from 3 to 6 (see e.g. Turner |T994| l. If that 
were the case for Westerlund 2, all the distance determinations 
in this paper would have to be revised downwards. To determine 
the value of Ry requires at least UBV photometry. Since we do 
not have our own new U -band measurements, we would have 
to rely on the data of Moffat et al. ( 119911 1 or Carraro & Munari 



(2004). Unfortunately, both data sets seem to have problems (see 
Sect. |2.1t and their U magnitudes disagree. Our attempts to use 
these archive photometric data failed to produce conclusive re- 
sults and we have to conclude that the reddening law cannot be 
quantified with the data currently available. 
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